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Abstract 

We study the phase diagram of a one-dimensional Hubbard model where, 
in addition to the standard nearest neighbor hopping t, we also include a 
next-to-nearest neighbor hopping t'. For strong enough on-site repulsion, this 
model has a transition at half filling from a magnetic insulator with gapless 
spin excitations at small t'/t to a dimerized insulator with a spin gap at 
larger t'/t. We show that upon doping this model exhibits quite interesting 
features, which include the presence of a metallic phase with a spin gap and 
dominant superconducting fluctuations, in spite of the repulsive interaction. 
More interestingly, we find that this superconducting phase can be reached 
upon hole doping the magnetic insulator. The connections between this model 
and the two chain models, recently object of intensive investigations, are also 
discussed. 
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I. INTRODUCTION 



The properties of correlated electrons confined to a double chain have recently attracted 
considerable attention both from the theoretical and experimental point of views. 

The theoretical analyses have been mainly focused on simplified models as two chains of 
electrons interacting via a short range repulsion (e.g. the Hubbard model) and coupled by 
a transverse hopping t±, or two t — J chains coupled both by a transverse hopping and by a 
transverse exchange J±. At half filling both models are equivalent to two coupled Heisenberg 
chains whose ground state has been found to be a spin liquid insulator with a gap in the 
excitation spectrum for arbitrary transverse coupling (!]-[§. Away from half-filling both 
models describe a metal which however maintains a finite gap for the spin excitations. This 
behavior suggests the existence of electron pairs which is confirmed by the evidence that the 
dominant fluctuations describe density waves and interchain-pairing fluctuations |2| §. 
The latter are expected to dominate for weak repulsion and sufficiently away from half filling 
or, in the t — J ladders, for strong J. 

From the experimental point of view, recent measurements on ladders compounds like 
SrCu2C>3 and (VO)2P20y [[U]] confirmed the theoretical prediction of a spin gap at half- 
filling. The transition upon doping from the spin-liquid insulator to the metal with a spin 
gap has also been verified experimentally in the Sr doped LaCu02.5 ]TT[] . Unfortunately no 
superconducting transition seems to occur down to 5K, which is however not in contrast 
with the theoretical predictions (it would imply either that the doping is still low or that 
the interaction is too strong). 

An important message which in our opinion arises from all the theoretical analyses of 
the two chain models and which is the subject of the present work, is that doping a ID spin 
liquid may indeed result in superconductivity also in the presence of repulsive interaction. 
The goal of this paper is to show that this feature is shared not only by two chain models 
but also by a wider class of ID models which do describe a spin liquid insulator at half 
filling. 
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Among the spin models which are known to exhibit a spin-gap in the excitation spectrum, 
a very simple and well studied model is the spin-1/2 Heisenberg chain with an additional 
next-to-nearest neighbor exchange 

L L 

Hjj> = J^Si- Si+i + J' ^2 Si - S i+2 . (1) 

i=l i=l 



If J' = this model is the well known Heisenberg model |f[2 |, which is characterized by 
gapless excitations and power-law decay of the spin correlations. If J' = J/2 the ground 
state is exactly known |13| and consists of a product of singlets among nearest neighboring 
sites (dimerized state). There are two of these states, which are related among each other 
by the translation of one lattice constant. A finite energy gap exists between these two 



degenerate states and the first excited ones fTM . The transition upon increasing J' from the 



gapless regime to the gapped dimerized state was studied using bosonization by Haldane [15 



who predicted the transition to occur at J' ~ J/6. Successively, Nomura and Okamoto fl6 
performed a detailed numerical investigation of the model and estimated a larger transition 
value of J' ~ J/4. 

A model of interacting electrons which in a particular limit reproduces the spin-model 
(H) is the Hubbard model with an additional next-to-nearest neighbor hopping (t — i! — U 
model), described by the Hamiltonian 

h = -tY.il ( c U +lCT + h.c) + f E E (4a- +2ct + H - c ) + ^X>T^i, ( 2 ) 

£T=T4 i=1 cr=T,l *=i »=i 

where Cj CT annihilates a spin a electron at site i and n io . = c\ a c ia . At half- filling and for 
U ^> (t, t') this model indeed maps onto (0) with J = At 2 /U and J' = At' 2 /U, and therefore 
it is a good candidate to study the properties upon doping of a spin liquid state. 

In this paper we study the phase diagram of (0) by making use of weak coupling Renor- 
malization Group (RG) and bosonization. We will show that, as a function of the parame- 
ters (electron density n, U/t and t'/t), the phase diagram is surprisingly rich. In particular 
we find that also in this simple case superconductivity may arise from doping the spin 
liquid insulator, even though the electron-electron interaction is repulsive. Moreover for 



t/4 < t' < t/2 we find a transition upon doping from a magnetic insulator at half filling to 
a metal with dominant spin and charge density wave fluctuations and finally to a supercon- 
ductor (for small U/t) or a metal with a spin gap and dominant dimer wave fluctuations (at 
larger U/t). Although the model is purely one dimensional, this behavior is quite suggestive 
especially for its similarity to the phase diagram of HTc superconductors. 



II. THE MODEL 

As stated in the Introduction, we are going to study the model described by the Hamil- 
tonian (0) which, in the absence of interaction, has the following energy dispersion relation 

e(k) = -2tcosk + 2t'cos2k. (3) 

Notice that the model has particle-hole symmetry if, at the same time, t' — > —t'. Let us first 
analyze the dispersion relation ([|) which is the starting point of our perturbative analysis. 

If if < i/4 the band minimum is at k = (see Fig.l). The model is then a simple one- 
band model and if the interaction U is turned on we expect a behavior qualitatively similar 
to the standard Hubbard model (if = 0). We are not going to discuss this case in much 
detail, since its behavior is very well known |T7 . 



If t' > t/4 the band minima ±fc m m move away from k = (which turns into a band 
maximum) and satisfy the relation 

COS rv m i n t . 

In this case two different situations may occur (see Fig.2): 

(1) if the density is such that the chemical potential is bigger than e(0) = — 2t + 2if, the 
model at low energy is effectively a one-band model, for which the previous conclusions 
for the case if < t/4 apply; 

(2) if, on the contrary, the chemical potential is smaller than e(0), there are four Fermi 
points (±kpi and ±/cp 2 ), thus the model at low energy behaves as a two-band model. 



At half filling this implies that: 

• if t' < t/2 there are only two Fermi points ±/ci? = ±7r/2. There is therefore a simple 
Umklapp scattering since Akp = 2n, exactly like in the standard Hubbard model; 

• if t' > t/2 there are four Fermi points (see Fig. 2) satisfying the relation 2kp 2 — 2k pi = 
tt. In this case, as we are going to discuss in the following Section, there is only a 
higher order Umklapp which involves four-electron scattering at the Fermi surface, 
since 4kp 2 — 4&fi = 2n. 

If U 7^ and one is interested in the low energy behavior, a standard approach for a ID 
system is to linearize the band around the Fermi points: e(k) = ±vp(k =F kp) (see Fig.3a) 
if there are only two Fermi points, while ei(k) = ^Vpi(k =p kpi) and e 2 (k) = ±vp 2 (k =F kp 2 ) 
if four Fermi points are involved (see Fig.3b). The linearization is assumed to be valid 
only within some cutoff range of width A. The interaction U causes scattering among these 
Fermi points. These scattering processes are relevant since they generate in perturbation 
theory logarithmic singularities (usual Cooper's singularity in the particle-particle channel 
and additional singularities in the particle-hole channels due to the nesting property of 
a ID Fermi surface). A standard way to cope with such a log-singularities is the weak 
coupling Renormalization Group (RG), which, along with the bosonization technique, is a 
very powerful tool in ID. Although our analysis will make use of these techniques, we are 
not going to introduce them since there exist a wide number of articles where they have been 
intensively discussed |L8|-fJ(J. For the two-band model, in particular, I will closely follow 



the analysis of Ref. 0. In this reference a two-band model resulting from a two-coupled 
chain model was analyzed both via bosonization and RG. The only difference with Ref. Q 
is that the inner bands (which is e\{k) in the present case and the anti-bonding band with 
transverse momentum k± = it in the two-chain model) have opposite slopes in the two cases. 
The correct mapping between the two models is therefore 0] 



On provision that the previous mapping is performed, the perturbation expansion of the 
t — t' — U two-band model and of the two-chain Hubbard model is exactly the same at low 
energy (apart from an important difference at half filling, see next Section). Therefore we 
can simply borrow all the results which have been obtained for the two-chain models and 



use them for the present case. This is what we are going to do in the following Sections pi 

Already at this stage it is apparent that the behavior of the two-chain models is similar 
to that of a single chain t — t' — U model, and that the feature which makes the two class 
of models equivalent is the presence of four Fermi points in some parameter range. 



III. THE MODEL AT HALF-FILLING 

If the density corresponds to one electron per site, two cases occur, as previously dis- 
cussed. 



A. t' < t/2 

If t' < t/2, the low energy model is a one band model with Fermi momenta ±7r/2. There 
is a relevant Umklapp which makes the system an insulator. However the spin excitations 
are gapless and, as a consequence, the spin correlations have a power law decay at large 
distance. The model, for what it concerns the spin degrees of freedom, behaves exactly like 
a Heisenberg model [T^.[T3 . 



B. t' > t/2 

If t! > t/2, the effective model is a two band model. It is therefore worthwhile to start 
with a broad outline of the behavior of such a two-band model in ID. 

Without Umklapp terms, two different phases exist depending on the ratio of the Fermi 
velocities Vpx/vpi and U/t. If t' ~ t/2, the Fermi velocity of the inner band vpi <C vp 2 . 
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In this case RG predicts 0,|2| that the model is a metal with four gapless excitations (two 
spin and two charge sound modes). The properties of the ground state can be inferred 
from the correlation functions which have the slowest decay at large distances. In this case 
these correlation functions describe spin and charge density waves at the incommensurate 
momenta 2kp 2 and 2kpi- By increasing t' also vpxjvpi increases and at a critical t' c a 



transition to a different phase occurs flj2lj22|1- In this new phase the model has a gap for 



the spin excitations, and a single gapless charge mode which corresponds to the ordinary 
zero sound @J^,[?j]. There are two competing correlation functions which have the slowest 
asymptotic decay. One is, in the two-chain language, the 4k f charge density wave |||],[7[]. 
In the language appropriate to the t — t' — U model this function translates into the Dimer 
Wave (DW) correlation function which decays at large distances like: 

. ^ , . ^ ... cos \2(kp2 — kpilx] cos(ttx) , jS 

Xdw(x) = (O DW (x)O DW (0)) 1 1 \ K ' J = — (4) 

the last equivalence being true only at half filling and 

DW {x) = S + {x)S-{x + a) - S + {x - a)S-(x) (5) 

being the dimer order parameter [ITH], with a the lattice constant. 



The other competing correlation function is what in the two-chain language has been 
identified as a kind of ci-wave superconducting correlation function (SC) [[| . In the t — t' — U 
model 

X sc(x) = (A(x)At(O)) ~ -L^, (6) 

where 

(7) 

P =± 

The Fermi operators ?/>'s in (|^) are defined around each Fermi point, i.e. 

ip P k F A x ) ~ eiphF ' X eikXc pk Fi +k,a, 

|fc|<A 



where i = 1,2. Notice that the existence of a spin gap already signals some kind of electron 
pairing. Due to the repulsive nature of the interaction, the pair wave function should have 
a minimum whenever the two electrons approach each other. This is accomplished by the 
minus sign in the expression of the pair operator A(x), Eq.(|7]), which in turns shows the 
importance of having more than two Fermi points at disposal. However, the existence of 
electron pairs does not necessarily imply dominant superconductivity. This depends upon 
the pair-pair interaction which in turn determines the value of the parameter K. 

From Eqs.(|)-(D, we see in fact that if K > 1/2 the pairing fluctuations indeed dominate 
over the DW fluctuations, while the opposite occurs if K < 1/2. According to bosonization 
T9| , |T2"|| , K is related to the charge compressibility. In particular 

1 4L d 2 E 



K ttv p dN 2 ' 



(8) 



where L is the length of the chain, E the ground state energy, N the electron number and 
v p the velocity of the charge zero sound. The latter can be determined numerically by 
calculating the energy gap between the ground state (for closed shells at total momentum 
P = 0) and the first excited state at total momentum P = 2n/L 



E{P = 2tt/L) - E(P = 0) = -^-v p . 

Lj 



The larger the electron-electron repulsion, the smaller the compressibility and consequently 
K, and the more unlikely is the dominance of superconductivity. 

Since we are at half filling, we have also to take into account Umklapp scattering. In 
this case there is only one higher order Umklapp process, which involves a four-electron 
scattering at the Fermi surface (see Fig.4). From dimensional analysis it turns out that this 
Umklapp is relevant if K < 1/2. In this case the zero sound mode acquires a gap and the 
model becomes insulating. Having discussed the possible phases of the t — t' — U model when 
the Fermi surface has two branches, let us study in detail their occurrence at half filling. 

If t' > t' c and K < 1/2, the Umklapp is relevant and therefore the model is insulating 
with a gap in the whole excitation spectrum and a finite average value of the dimer order 
parameter Eq.@ 
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(O dw {x)) = (-l) x • const. 

This insulating phase certainly occurs for U ^> (t, t') when the mapping to the spin model 
(|I|) is justified. On the other hand, for very small U, the parameter K can be evaluated 
by perturbation theory and it turns out to be close to one, modulo corrections of order 
U. Therefore, provided perturbation theory is valid, K > 1/2 for U <C t, which implies 
that the Umklapp scattering is irrelevant and the model is metallic with the dominant 
superconducting fluctuations Eq.(H). Consequently we expect a transition at a finite U = U c 
from a metal with superconducting correlations directly to an insulator with a dimer order. 
In Fig. 6 we have drawn a qualitative phase diagram for if > t' c > t/2. At half filling, n — 1 
in the figure, there is a critical U which separates the insulating regime (the bold line in 
the figure which we label DI, meaning a dimer insulator) at larger U from the metal with 
superconducting fluctuations at smaller U (which we label SC). 

If t/2 < t' < t' c , we still expect a metal-to- insulator transition at a finite U, but this time 
the metal has no spin gap and shows dominant density wave fluctuations. The properties 
of the insulating phase into which the above metal transforms at large U can not be simply 
deduced by means of RG, whose validity is doubtful at finite U. However, we tend to believe 
|23| that this insulator should have the same properties of the dimer insulator which occurs 
for t' > t' c . 

Notice that the behavior of the t — t' — U model at half filling is different from the 
behavior of the two chain models also at half filling. There, the Umklapp term is a two- 
electron scattering process and is relevant for any K < 1, which implies that the model is 
an insulator for any f/^ @. 

To conclude this section, we like to point out that, according to our analy- 
sis, the transition at large U between the insulating phase with power lay de- 
cay of the spin correlations and the dimer insulator with a spin gap is predicted 
to occur at t' ~ t/2 or, in terms of the exchange couplings, at J' ~ J/4. 
This is exactly the value found by numerical investigation of the spin model (JI|) in Ref. [|16|| . 
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This coincidence, which might well be accidental, is quite surprising, since our prediction is 
based simply on band-structure arguments (modification of the Fermi surface). In fact, one 
would rather believe that band structure details are irrelevant for electron systems in the 
strong correlation limit where the interaction is much larger than the bandwidth, which is 
the case of the model ([!]). 

IV. THE MODEL AWAY FROM HALF FILLING 

We have seen that already at half filling the t — t' — U model shows the unusual property 
of a transition at a critical U from a metal with superconducting fluctuations to an insulator, 
when t' > t/2. Away from half filling the behavior is even more interesting, for smaller and 
larger t'. Let us consider in detail the various possible scenarios. 

A. t' < t/4 

In this case the model is for any filling always an effective one-band model, where nothing 
special occurs. Whether we dope with holes or electrons, as soon as we move away from half- 
filling the charge gap closes and the system becomes a metal with gapless spin and charge 
sound modes. The dominant fluctuations are both spin and charge density wave fluctuations 
(SDW and CDW) at momentum 2kp . In Fig. 5 we have drawn the phase diagram for this 
case as a function of U/t and of the density n. The bold line at density n — 1, labeled MI, 
identifies the magnetic insulator with power law decay of the spin correlations, while the 
rest of the phase diagram has been labeled with SDW/CDW implying it is a metal with 
dominant density wave fluctuations. These labels will be used with the same meaning also 
in the following cases. 
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B. t' > t/2 



Here, it makes a difference whether we dope with holes or electrons (the model is not 
particle- hole symmetric). In Fig. 6 we have drawn a qualitative phase diagram for t' > t' c > 
t/2. 

For hole doping (n < 1), the effective low energy model involves always two bands. 
Therefore we predict that for any hole doping the spin gap will survive. As concerns the 
charge gap, it will immediately disappear as soon as we move away from half filling. There- 
fore, as U increases, we expect a crossover from a metal with dominant superconducting 
fluctuations [see Eq.@], which we still label in Fig. 6 as SC since it is continuously con- 
nected to the analogous state at n — 1, to a metal with dominant dimer wave fluctuations 
[see Eq.flD], which is labeled by DW. 

At this point, it is worthwhile to discuss briefly the properties of the weakly doped 
dimer insulator. If we approach half-filling with U > U c then, according to the theory 
of the incommensurate-to-commensurate transitions in ID p4[ , K tends asymptotically to 
the value 1/4 and, exactly at half-filling, it jumps abruptly to zero. With the asymptotic 
value K = 1/4 valid at low doping, the dimer wave correlation function Eq.(^) decays 
like Xdw ~ ( — 1) X /V^) while the superconducting correlation Eq.(|6|) decays quadratically 
Xsc ~ 1/x 2 . Notice that these power law decays are typical of the Green functions and of 
the density-density correlation functions, respectively, of a hard core Bose gas. The situation 
is opposite for the metallic phase at U <C t. In this case, provided perturbation theory is 
valid, K ~ 1 and the behavior of the two correlation functions are exchanged: xsc ~ V\/^ 
while xdw ~ {—l) x /x 2 . 

In the case of electron doping the situation is in general different but for low doping 
where all is the same except the chemical potential which moves up instead of down (see 
Fig. 2). At the same time the Fermi velocities vfi and vf2 of the two linear bands get more 
and more different (actually vfi — > 0). As we said in the previous Section, at a critical 
value of vf 2 /vfi, or equivalently a critical density, n cl in Fig.6, RG predicts a transition to 
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another phase where also the spin gap closes (see e.g. Appendix B in Ref. (4j). In this phase 
the model is a metal with four gapless sound modes (two spin and two charge modes). The 
dominant fluctuations describe charge and spin density waves. We have labeled this phase 
in Fig. 6 as SDW/CDW II, implying that the number of gapless excitations is twice that of 
the phase SDW/CDW. 

Finally, at a second critical doping n C 2, the topology of the Fermi surface changes from 
a four-point to a two-point Fermi surface. At the transition the density of states of the 
inner band diverges due to a van Hove singularity. For this reason we are not able to 
predict what happens exactly at the transition. According to Balents and Fisher ||, the 
van Hove singularity induces again a spin gap and therefore they expect the properties of 
the model to be similar to those at low doping. On the other hand, if we assume that 
the RG equations of Ref. || can be extended up to a very large Vp 2 /vpi (where their 
validity is not fully guaranteed), we would rather expect that for Vf^I^fi 3> 1 the two linear 
bands effectively decouple. In this case the transition would be a standard metal-to-metal 
transition with a topological modification of the Fermi surface. The van Hove singularity 
related to the low (hole) doping of the inner band is not expected to play any fundamental 
role, similarly to what happens to any one-band model close to filling zero or one. Although 
we have no rigorous proof, we tend to believe that nothing special occurs at the transition 
(the latter scenario), rather than in the scenario proposed in Ref. 0. Coming back to the 
phase diagram, in the region where only two Fermi points are involved the model should 
be metallic with still dominant density wave fluctuations (see Fig. 6), and only two gapless 
modes (therefore this phase is labeled SDW/CDW). 

In Fig. 6 we have assumed that the critical density n c \ decreases by increasing U. This is 
true for very small U where perturbation theory is valid. For finite U, unfortunately, we have 
no reliable method to evaluate n c \. It might well be that n cl — ► n c2 at U (t,f), which 
would be more consistent with the proposed transition upon increasing U at half filling and 
for t' ~ t/2 from the phase CDW/SDW II to the dimer insulator DI. The uncertainty about 
the precise behavior of function of U is also the reason why we have prefered not 
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to draw even a qualitative phase diagram for t' close to t/2. 

C. t/A < t' < t/2 

This case is in our opinion the most interesting one for its surprising similarities with the 
behavior of HTc compounds. We said in the previous Section that at half filling the model 
is an insulator (MI in Fig.7 at n = 1) with power law decaying spin correlations. If we dope 
with electrons (n > 1), the charge gap suddenly closes and the model turns into a metal 
with dominant spin and charge density wave fluctuations, for arbitrary doping (SDW/CDW 
in Fig.7). 

If we dope with holes (n < 1), at low doping the Fermi surface still consists of two points, 
and therefore the properties are similar to what we just described for electron doping, i.e. 
the phase is a SDW/CDW. However, as the chemical potential moves down (see Fig. 2), the 
same situation encountered earlier re-appears, now in reverse. At a lower critical doping 
(n c2 , in analogy with the previous case) there will be a topological modification of the Fermi 
surface from a two-point to a four-point surface. The model turns from a one-band to a two- 
band model. Until the Fermi velocities remain quite different (vpz/vpi ^> 1) the system is a 
metal with four gapless modes (two charge and two spin modes) and dominant density wave 
fluctuations (phase SDW/CDW II in Fig.7). At a critical value of Vp 2 /vpi, or equivalently 
a critical doping n cl , a spin gap opens and only one charge mode (the ordinary zero sound) 
remains gapless. The dominant fluctuations are dimer waves Eq. (|5]) or pairing fluctuations, 
Eq.(^). The one which dominates depends on the value of the parameter K Eq.(|8|), which in 
turns depends on U. As we said before, if K < 1/2 the dimer fluctuations win (DW phase 
in Fig.7), while for K > 1/2 superconductivity is more relevant (SC phase in Fig.7). Like 
in the previous case, we can not establish the precise behavior at large U of the critical line 
between the CDW/SDW II phase and the SC or DW phases. Therefore the shape of the 
phase diagram drawn in Fig.7 has not to be taken too literally close to that critical line. 

Notice that the appearance of the spin gap upon hole doping the magnetic insulator 
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indicates that in this model the holes increase the spin frustration. In the large U limit, one 
can map the t — if — U model onto a generalized t — J model with the Hamiltonian 

# = E E (cU+i CT + H.c.) + f E E (cU +2ff + H.c.) 

L L 
+ J E "^i ' + J' E ' *S''i+2, 

defined in the reduced Hilbert space where double occupancies are forbidden. This model 



with t' = has been numerically investigated by Ogata, Luchini and Rice ||25|| . They find 
that the hole doping effectively reduces the frustration due to J'. Since we instead find 
an increase of frustration, we have to conclude that this is mainly a consequence of the 
next-to-nearest-neighbor hopping t' [2B|. 



V. CONCLUSIONS 

In this paper we have studied the phase diagram of a one-dimensional t — t' — U model 
where, in addition to the on site repulsion U and nearest neighbor hopping t, we have 
included a next-to-nearest neighbor hopping t'. Although very simple, this single chain 
model has the interesting property to be at half filling and large U either an insulator 
with gapless spin excitations or a dimerized insulator with a spin gap ||15|| . The transition 
between the two insulators should occur, according to our analysis, at t' ~ t/2. This model 
is therefore particularly suited to study the occurrence of superconductivity upon doping an 
insulator which has a charge gap and a spin gap. In fact, recent theoretical and numerical 
investigations of two coupled chain models PH, suggest that the presence of a spin gap in 
the insulating phase of two chains at half filling may lead to superconductivity upon doping 
and for not too large repulsion. We have shown that precisely this behavior is realized 
here. Moreover, the weak coupling analysis seems to suggests that the key feature which is 
responsible for the spin gap and possibly for the superconductivity is the topology of the 
Fermi surface in both models, which in the interesting parameter range has two branches, 
i.e. four Fermi points. 
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The present t — t' — U model has additional properties which make this model interesting 
in its own right. In particular at half-filling and for t' sufficiently larger than t/2, we predict 
a direct transition at a critical U c from a metal with dominant superconducting fluctuations 
at U < U c to a dimerized insulator at U > U c . 

Interestingly, for i/4 < t' < t/2 the phase diagram for hole doping shows some similarities 
with the phase diagram of HTc compounds (see Fig. 7). At half-filling we have a magnetic 
insulator with power law decay of the spin correlations. For low doping we move to a metal 
with dominant spin and charge density waves fluctuations. Above a critical doping, a spin 
gap opens and the model has either dominant superconductivity or dimer waves depending 
upon the strength of the on site repulsion. 
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FIG. 2. Energy dispersion relation of the t — t' — U model for t' > t/4. Also drawn are the 
chemical potentials corresponding to two different fillings: refers to the case when only one 

(2) 

band is involved at low energy while e F refers to the case when two bands are involved. 
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FIG. 3. The effective low energy models: (a) simple one-band model; (b) two-band model. 




FIG. 4. The relevant Umklapp scattering when the Fermi surface is made by four Fermi points. 
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FIG. 5. Phase diagram of the t — t' — U model for t' < t/4 as a function of f//t and of the 
density n. 
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FIG. 6. Qualitative phase diagram of the t — t' — U model for t' sufficiently larger than t/2 as 
a function of U/t and of the density n. 
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FIG. 7. Qualitative phase diagram of the t — t' — U model for t/4 < t' < t/2 as a function of 
U/t and of the density n. 
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